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UPPER BOUNDS FOR THE ATTRACTOR 
DIMENSION OF DAMPED NAVIER-STOKES 
EQUATIONS IN 

ALEXEI ILYINi, KAVITA PATNI ^ AND SERGEY ZELIK^’^ 

Abstract. We consider finite energy solutions for the damped 
and driven two-dimensional Navier-Stokes equations in the plane 
and show that the corresponding dynamical system possesses a 
global attractor. We obtain upper bounds for its fractal dimension 
when the forcing term belongs to the whole scale of homogeneous 
Sobolev spaces from —1 to 1. 


1. Introduction 

The theory of global attractors for the 2-D Navier-Stokes system 

J dtU + (m, Va;)M + + g, 

“O’ divM = 0. 

has been a starting point of the theory of inhnite dimensional dis¬ 
sipative dynamical systems and remains in the focus of this theory, 
see n [ini HI na EEi EH EE] and the references therein. 

In the case of a bounded domain the corresponding dynamical 
system has a global attractor in the appropriate phase space. The 
attractor has hnite fractal dimension, measured in terms of the dimen¬ 
sionless number G (the Grashof number), G = 

The best known estimate in the case of the Dirichlet boundary con¬ 
ditions M|ao = 0 is (see 1371 ) 

dimjA^CDG, (1.1) 

while in the case of a periodic domain x G [0,27rL]^ the estimate can 
be signihcantly improved (see my 

dim/A < CperG^'^^(ln(l -b G))^/^, 

and, moreover, this estimate sharp up to a logarithmic correction as 
shown in [2H] , see also HT] for the alternative proof of the upper bound. 
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Finding explicit majorants for the constants cd and Cper amounts to 
finding sharp or explicit constants in certain Sobolev inequalities and 
spectral Lieb-Thirring inequalities. For example, cd < [23] . 

and the majorant for Cper can be easily written down using the recent 
result [7] on the sharp constant in the logarithmic Brezis-Gallouet in¬ 
equality (which is essential for the attractor dimension estimate in the 
periodic case). 

In unbounded channel-like domains the Navier-Stokes system with 
Dirichlet boundary conditions is still dissipative in view of the Poincare 
inequality. In particular, if the case of hnite energy solutions is con¬ 
sidered, the associated semigroup possesses a compact global attractor 
similarly to the case of bounded domains. Up to the moment, there are 
two alternative ways to establish this fact. The hrst one is based on 
the weighted energy estimates and careful analysis of the Leray projec¬ 
tion in weighted Sobolev spaces, see [21 [29] and references therein and 
the second one utilizes the so-called energy method (which will be also 
used in our paper) and the energy equality, see [3, |Ml EZ] • 

However, in contrast to the case of bounded domains, the solution 
semigroup is no longer compact, but is rather asymptotically compact 
and this fact strongly affects the existing upper bounds for the dimen¬ 
sion of the attractor. Indeed, the best known estimate for the case of 
channel-like domains obtained in [3l| can be written as follows 


dim/ A < 


2 Xlu^' 


Here Ai > 0 is the bottom of the spectrum of the Stokes operator in 
the channel, and clt is the universal Lieb-Thirring constant, see fl3.3p . 
For example, in a straight channel of width d we have Ai > so 

that in this case we obtain 


dim/ A < 



d‘hf 

,.4 


We observe that these estimates are proportional to (unlike fll.ip 
which is proportional to On the other hand, it is worth men¬ 

tioning that, to the best of our knowledge, no growing as z/ —)■ 0 lower 
bounds for the dimension of the attractor are known for the case of 
Dirichlet boundary conditions, regardless whether the underlying do¬ 
main is bounded or unbounded. So, in contrast to the case of periodic 
boundary conditions, the behaviour of the attractor’s dimension as 
f ^ 0 remains unclear for the case of Dirichlet boundary conditions 
even in the case of bounded domains. 
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We mention also that keeping in mind the Poiseulle flows and the 
mean flux integral, it seems more natural to consider the infinite energy 
solutions for the Navier-Stokes system in a pipe. In this case, the 
system remains dissipative and the existence of the so-called locally 
compact attractor can be established, see Elia. The dimension of this 
attractor may be infinite in general, but will be finite if the external 
forces g{x) decay to zero as |a:| —>■ oo (no matter how slow this decay 
is), see [221112] and the references therein. 

In the whole the Laplacian is not positive-definite and the Navier- 
Stokes system is not dissipative at least in a usual sense even in the case 
of zero external forces and finite energy solutions, see e.g. [35]) see also 
[211130] and the references therein concerning the decay properties of 
various types of solutions for the Navier-Stokes problem with zero ex¬ 
ternal forces. The presence of external forces makes the problem more 
complicated and usually only growing in time bounds for the solutions 
are available. We mention here only the recent results concerning the 
polynomial growth in time for the so-called uniformly local norms of 
infinite-energy solutions obtained in [221 113], see also the references 
therein. 

Let us now consider the damped and driven Navier-Stokes system 

I dtu -F {u, V^)u V^p + au = lyA^u + g, 
divM = 0. 

with additional dissipative term au. The drag/friction term au, where 
a > 0 is the Rayleigh or Ekman friction coefficient (or the Ekman 
pumping/dissipation constant), models the bottom friction in two- 
dimensional oceanic models and is the main energy sink in large scale 
atmospheric models [3T] . 

The analytic properties of system fll.2l) (such as existence of solu¬ 
tions, their uniqueness and regularity, etc.) remain very close to the 
analogous properties of the classical Navier-Stokes equations or Euler 
equations if the inviscid case i/ = 0 is considered. However, the fric¬ 
tion term au is very essential for the dynamics since it removes the 
energy which piles up at the large scales and from the mathematical 
point of view compensates the lack of the Poincare inequality. This 
makes the Navier-Stokes system and even the limit Euler system dissi¬ 
pative whatever the domain is and allows to study its global attractors 
in various phase spaces. For instance, the so-called weak global at¬ 
tractor for the inviscid case is constructed in [22] for the case of finite 
energy solutions; its compactness in the strong topology related with 
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the if^-norm is verified in m and na for the cases of finite and infi¬ 
nite energy solutions respectively; the inviscid limit z/ —)■ 0 is studied 
in |T6] including the absence of the so-called anomalous dissipation of 
enstrophy; the existence of a locally compact global attractor in the 
uniformly local phase spaces is established for the viscous case u > 0 
in [13|; see also 0. ini for the existence and uniqueness results for the 
stationary problem and the stability of stationary solutions for fll.2l) 
with z/ = 0. 

From the point of view of the attractors and their dimension the sys¬ 
tem fll.21) in the case of the periodic domain a: G [0, 27rL]^ was studied 
in [21], where it was shown that the corresponding dynamical system 
possesses a global attractor A (in whose fractal dimension is finite 
and satisfies the following estimate 

(1,3) 

\ ua 8 J 

where the values of the constants are updated in accordance with [2S] . 
The first estimate is enforced in the regime v/L?‘ 3> a, while the second 
estimate is enforced in the opposite regime z//L^ -C a. We observe 
that both estimates are of the order 1/z/ as z/ —)■ O’*" if all the remaining 
parameters are fixed. It was also shown in [2l| that this rate of growth 
of the dimension is sharp, and the upper bounds were supplemented 
with a lower bound of the order 1/z/, based on the instability analysis of 
generalized Kolmogorov flows. The finite-dimensionality of the global 
attractor in the uniformly local phase spaces under the assumption 
that the external forces g{x) decay to zero as |a:| —)• cx) has been proved 
recently in [32|, but no explicit upper bounds for the dimension was 
given there. 

We would also like to point out that starting from the paper [27] the 
Lieb-Thirring inequalities are an essential analytical tool in the esti¬ 
mates of global Lyapynov exponents for the Navier-Stokes equations. 
This fully applies to our case. 

We now observe that the first estimate in fll.3l) blows up as the 
size of the periodic domain L —)■ cxo. On the other hand, the second 
estimate survives (the homogeneous if^-norm is scale invariant in two 
dimensions). Therefore, one might expect that this estimate holds for 
L = oo, that is, for a; G and a motivation of the present work is to 
show that this is indeed the case. 

In this paper we study the damped and driven Navier-Stokes system 
fll.211 in in the class of finite energy solutions and our main aim is 
to obtain explicit upper bounds for the attractor’s dimension in terms 
of the parameters z/ and a and various norms of the external forces g. 
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In section [2] we recall for the reader convenience the proof of the 
well-posedness and derive the energy equality. Then using the energy 
equality method [301131] we establish the asymptotic compactness of the 
solution semigroup and, hence, the existence of the global attractor A. 

In section [3] we consider the case when the right-hand side g belongs 
to the scale of homogeneous Sobolev spaces H^, s G [—1,1] and derive 
the following estimate(s) for the fractal dimension of the attractor A: 


dim f A < -^ 

^ “ 64\/3 


/ i+igi Y^i 1 

\1 — |s|/ 



■S ^ [“!) !]• 


In particular, for s = 1 we obtain 


dim f A < 


1 II curl 
16x/3 


which up to a constant coincides with the second estimate in fll.3p . 
proving thereby our expectation. 

In this paper we use standard notation. The L 2 -norm and the cor¬ 
responding scalar product are denoted by || ■ || and (•,•)• 


2. A PRIORI ESTIMATES, WELL-POSEDNESS AND ASYMPTOTIC 

COMPACTNESS 

We study the damped and driven Navier-Stokes system fll.2p in 
Here u{t, x) = {u^, v?) is the unknown velocity vector held, p is the un¬ 
known pressure, g{x) = {g^, g'^) is the given external force (and without 
loss of generality we can and shall assume that div = 0), the advection 
term is 

2 

(u, Vx)v = 

i=l 

and a > 0, u > 0 are given parameters. We restrict ourselves to 
considering only finite energy solutions, so we assume that 

g G [L^(R^)]^ Uoe'H:= {uq G divuo = 0}, 

and by dehnition u = u{t, x) is a weak solution of fll.21) if 

M G C'(0,T;?f) nL2(0,T; [Tr^(M2)]2)^ y > 0, (2.1) 


and fll.2p is satished in the sense of distributions. This means that 
for every divergence free test function (p(t,x) = div(p = 0, 
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G C^{R+ X M^), the following integral identity holds 


— / {u,dt(j))dt+ / {{u^Vx)u, (f) dt + a {u,(p)dt+ 


+ u {VxU,Vx^)dt = / {g,^p)dt. (2.2) 


The following fact concerning the global well-posedness of the Navier- 
Stokes equations in 2D is well-known, see, for instance. 


Theorem 2.1. For any uq G R there exists a unique solution u of prob¬ 
lem fll.2p satisfying fl2.ip and fl2.2p . Moreover, this solution satisfies 
the following energy equality for almost all t > 0: 


1 d 

2 dt 




(s.'“(*))■ 


(2.3) 


Proof. For the convenience of the reader, we remind below the key 
steps of the proof. It is strongly based on the so-called Ladyzhenskaya 
interpolation inequality 

II’^IIl 4(R2) ^ ^ll'^ll ll^a:Il|| (2-4) 

which holds for any u G Indeed, this inequality together with 

fl2.ip implies that any weak solution u satishes 

u G L^(0,T; 


and, integrating by parts and using the Cauchy-Schwartz inequality, 
we see that 

2 

|(m, Vx)m,(P)| = \ < C\\u\\l4^\\m. 

i,j=^ 


Therefore, 


/ {{u,VxU),gFj dt < C||-u||^4(0,T;L4(]r2))||95||l2(0,T;H1(R2)) 

Jr 

and 

{u,Vx)u G = [L\t),T-V})]\ (2.5) 

where, as usual, R} is a subspace of [iL^(R^)]^ which consists of diver¬ 
gence free vector helds and R~^ := \R}]*. Thus, from fl2.2p . we see 
that dtu G L^(0,T;'H“^) as well, and all terms in fl2.2p make sense for 
the test function (p = m G L^(0, T; "H^), so by approximation arguments 
(and the fact that the divergent free vector helds with compact sup¬ 
port are dense in R) we may take (p = m in fl2.2l) . Then, using the 
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well-known orthogonality relation 

((m, Vcc)v, v) = -^(div u, = 0, (2.6) 

we end up with the desired energy equality (12.31) . see, for instance, [37] 
for the details. 

The existence of a weak solution can be obtained in a standard way 
based on the energy equality (12.3p either directly by the Faedo-Galerkin 
method or by approximation the problem (11.21) in by the correspond¬ 
ing problems in bounded domains (see e.g., [271EB] for the details), so 
we leave this proof to the reader and remind below the proof of the 
uniqueness. Let ui and U 2 be two weak solutions of problem (11.21) and 
let V = ui — U2. Then v solves 


dtv -t- (ui, V 2 ;)n -F (n, Vx)u2 + av + VxV = v/S.xV, div n = 0. 


Multiplying this equation by v and integrating over a: G , analogously 
to the energy equality, we have 


1 d 

2 dt 




-((«!, Vjtl), tl) - ((tl, VJtii, k). 


(2.7) 

The hrst term on the right-hand side of this equality vanishes in view 
of (12.61) and the second term can be estimated using the Ladyzhenskaya 
inequality as follows: 


|((n, Va;)M 2 ,p)| < ||n||^4||Va;M2|| < G||n||||V,j,n||||Va:M 2 || < 

< Z/||Va;T||^ + Va;tl2|n|T||^- 

Inserting this estimate into the right-hand side of (12.71) we obtain 

< C^iy-^\\VxU2{t)f\\v{t)f, 
and the Gronwall inequality gives 

\Mt)-U2it)f < -n2(0)|p. 

Thus, the uniqueness is proved and the theorem is also proved. □ 


Corollary 2.2. The weak solution u{t) of problem (II.2p satisfies the 
following dissipative estimate: 


2 v 



\u{t)r< 


-f 0-2 


Vxu{s)fds<a-^T\\gr + \\u{t)f- 


( 2 . 8 ) 
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Proof. Using the Cauchy-Schwartz inequality, the hrst estimate follows 
from 02.31) by the Gronwall inequality, and the second estimate is proved 
by integrating 02.3p . □ 

Thus, equation 01.21) dehnes a solution semigroup S'(t) ; "H -A "H; 

S(t)uo := u(t), Uq eH 

where u{t) is a weak solution of equation 01.2p with the initial data 
m(0) = Uq. Moreover, this semigroup is dissipative according to esti¬ 
mate 02.81) and is Lipschitz continuous: for any Mi, M 2 G 

||S'(t)Mi - S{t)u2\\ < Ge^*||Mi - M 2 II, (2.9) 

where the constants C and K depend only on ||mq||. Our next aim is 
to verify the existence of a global attractor for this semigroup. For the 
reader convenience, we hrst remind its dehnition, see [U HOl HH [37] for 
more details. 

Definition 2.3. Let S{t), t > 0, he a semigroup acting in a Banach 
space H. Then the set M C "H is a global attractor of the semigroup 
S{t) if 

1) The set A is compact in "H. 

2) It is strictly invariant: S{t)A = A. 

3) It attracts the images of bounded sets in as t —)■ 00 , i.e., for 
every bounded set B gH and every neighborhood 0{A) of the set A 
in PL there exists T = T{B,0) such that 

S{t)B c 0{A) 

for all t >T. 

To state the abstract theorem on the existence of a global attractor, 
we need more dehnitions. 

Definition 2.4. Let S{t) be a semigroup in a Banach space Pi. Then, 
a set B C Tf is an absorbing set of S(t) if, for every bounded set B G PL, 
there exists T = T{B) such that 

S{t)B G B. 

A semigroup S{t) is asymptotically compacts for any bounded sequence 
Mg in B and for any sequence tn —t 00 , the sequence S{tn)uQ is pre¬ 
compact in B. 

In order to verify the existence of a global attractor we will use the 
following criterion, see [H, [10], [26] [36] EZ] for its proof. 
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Proposition 2.5. Let S{t) be a semigroup in a Banach space H. Sup¬ 
pose that 

1) S{t) possesses a hounded closed absorbing set B C H; 

2) S{t) is asymptotically compact; 

3) For every fixed t > 0 the map S(t) : B ^ Ti is continuous. 

Then the semigroup S(t) possesses a global attractor A C B. More¬ 
over, the attractor A has the following structure: 

where K. C is the set of complete trajectories m : R —)■ "H o/ 

semigroup S{f) which are defined for all t eM. and hounded. 

The next theorem which establishes the existence of a global attrac¬ 
tor for the solntion semigronp S{t) associated with eqnation (ll.2|l is 
the main resnlt of this section. 

Theorem 2.6. The solution semigroup S{t) of the damped Navier- 
Stokes problem fll.2p possesses a global attractor A in FL. 

Proof. We will check the assnmptions of Proposition 12.51 Indeed, the 
hrst assnmption is satished due to the dissipative estimate f|2.8p and 
the desired bounded and closed absorbing set can be taken as 




( 2 . 10 ) 


The third assumption is also satished due to estimate fl2.9p . Thus, 
we only need to check the asymptotic compactness. We will use the 
so-called energy method (see O [30l [M]) in order to do this. 

Indeed, let G "H be a bounded sequence of the initial data. Then, 
due to estimate fl2.8p . we may assume without loss of generality that 
E B. Let Unit), t > —tn, tn -E oo, be the sequence of solutions of 

dt^n T x)^n T OCUji -|- ^xPn T g, 



diVM„ = 0, Un 


Then, Mn(0) = S{tn)u^ and we only need to verify that the sequence 
{m„(0)}5^o is precompact in TL. In order to do so, we hrst verify that 
there exists a subsequence (which we also denote by Un for simplicity) 
such that 


M„(0) ^ m(0) 


( 2 . 11 ) 


converges weakly in TL to m( 0), where u{t) is a complete bounded tra¬ 
jectory u E 1C. We hrst note that due to the dissipative estimate fl2.8p . 


Un\\L^{-T,T-,'H) + SUp \\Un\\L'^(t,t+l-,H^) ^ C, 


t>-T 


( 2 . 12 ) 
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where T < tn and C is independent of n and T. Moreover, from the 
Ladyzhenskaya ineqnality f|2.4p we conclnde also that Un is bonnded in 
L‘^{—T, T ; L^). Thus, by the Banach-Alaoglu theorem, we may assume 
without loss of generality that 

u"" ^ u weak-star in T ; "H) 

and weakly in L^{-T, T; H^) n L^(-T, T; L^) 

for every T G N to some function u G L°°(R; H^) which also 

satishes estimate fl2.12p . Moreover, analogously to fl2.5p . we conclude 
that dtUn is uniformly bounded in t + 1; Thus, without loss 

of generality 

dtu^ dtu 

in the space t -|- 1; 'H~^) for all f G M. Using the embedding 

L^{t, t + 1; n~^) n L^{t, t + 1;c c{t, t + i-,n), 

see m Lemma IIL1.2], we see that Mn(0) ^ m( 0). Thus, to verify 
fl2.1ip . we only need to check that u{t), t G M is a weak solution of 
fll.2D . In other words, we need to pass to the limit as n —)■ oo in the 
analogue of fl2.2p : 

- / (Un, dt(j)) dt + / {{Un,Vx)Un,(p)dt + a / {Un, (f) dt+ 

t/ M M 

+ u f {VxUn,'Vx(p)dt = f {g,(p)dt, (2.13) 

M t/ M 

where ^p{t,x) G is an arbitrary fixed divergence free function. 

Since passing to the limit in the linear terms is evident, we only need to 
pass to the limit n —)■ oo in the non-linear term. Integrating by parts 
and rewriting the nonlinear term in the form 

2 

*,i=i 

and using the fact that the support of ip is hnite, we conclude that, for 
passing to the limit in the nonlinear term, it is enough to verify that, 
for every hxed R,T > 0, 

Un ^ u strongly in L^(—T, T; L^(S(f)), (2.14) 

where Bq stands for the ball of radius R in centered at the ori¬ 
gin. To verify the convergence fl2.14p . we recall that the sequence Un 
is bounded in the space L‘^{—T,T-,R^{Bq)) due to the dissipative esti¬ 
mate fl2.8p . Moreover, analogously to fl2.5p but using the test functions 
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we see that dtUn is bounded in the space 
L‘^{-T,T-H-\B^)). Thus, since 

h'(bS)) c h{b^) c «-■«) 

and the hrst embedding is compact, the compactness theorem (see, for 
instance [33, Theorem III.2.1]) implies that the embedding 

is compact. Therefore, is precompact in L‘^{—T, T ; L‘^{Bq)) for every 
i? > 0, T > 0 and passing to a subsequence if necessary, we conclude 
that the convergence fl2.14l) indeed holds. Thus, passing to the limit in 
the nonlinear term of fl2.13p is verihed and u is a weak solution of fll.2l) 
which is dehned for alH G M and bounded, so u G /C. This means that 
the convergence 02.111) is verihed. 

We are now ready to verify that 

Mn(0) ^ m( 0) strongly in "H (2-15) 

and hnish the proof of the theorem. We multiply the energy equal¬ 
ity 02.31) for the solutions Un by and integrate from —tn to 0: 


|u„(0)||^ = -2 / e^'^^\\V,Un{s)rds+ 

J —tn 


+ ||Mn(-tn)H e 2 / (2.16) 


tn 


We want to pass to the limit n —)■ cx) in this equality. Indeed, using 
the weak convergence —)■ u in H^) implying that 


lim sup —2 




tn 


= —2 lim inf 

n—>-cxD 


e"“1|V^Un(s)|rds < -2 / e"“"||V^u(s)||"ds 


— tn 


and the uniform bounds 02.121) . we see from 02.161) that 


limsup ||Mn(0)||^ < —2 


^2as I 


V 2 .M(s)|r ds + 2 


e '^"ig,u{s))ds. 


On the other hand, thanks to the energy equality, for the whole-line 
L^-bounded solution u G /C we have 


■*^(0)11 = “2 / Va;u(s)|f ds-I-2 / {g, u{s)) ds 
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and, therefore, taking into the acconnt the weak convergence (12.lip , 
we hnally arrive at 

limsnp ||Mn(0)|p < ||■u(0)||^ < liminf ||■u„(0)|p. 

n—>-oo n—>-oo 

Thns, lim^^oo ||'Wn(0)|| = ||m( 0)||, and the strong convergence fl2.15p is 
proved: 

lim ||n„(0 )-m( 0)|P = lim (||m„(0)||^ + ||m(0)||^ - 2(m„(0),m(0))) = 0, 

n—>-oo n—>-oo 

and the theorem is also proved. □ 


To conclnde the section, we also discnss the extra regnlarity of the 
attractor A. We say that n : R —)■ ^ is a strong solntion of the damped 
Navier-Stokes eqnations fll.2p if 

u e ^(0,T; n^) n L2(0,T; ^^(R^)) 

and eqnation fll.2p is satished in the sense of distribntions. The follow¬ 
ing analogne of Theorem 12.11 gives the global well-posedness of strong 
solntions of the Navier-Stokes problem. 


Theorem 2.7. Let Uq G . Then, the weak solution u constructed in 
Theorem \2.1\ is actually a strong solution and satisfies for almost all 
t > 0 the following analogue of the energy equality: 


1 d 

2 dt 


VxM(t)||^ + p||Aa,M(t)||^ -F a\\V xU{t)\\‘^ 


{g,Axu{t)). 


(2.17) 


The proof of this theorem is analogous to Theorem 12.11 and even sim¬ 
pler since the solution u is a priori more regular now, see, for instance, 
[37] for more details. We only mention here that the identity fl2.17|) fol¬ 
lows by multiplication of equation fll.2p by A^u, integration over a; G R^ 
and using the well-known orthogonality relation {{u,'Vx)u, A^u) = 0 
which holds for any u G [if^(R^)]^, divu = 0. In fact, integrating by 
parts and setting oj := curln = dx^u^ — dx 2 U^, the divergence theorem 
gives 

((u, 'Vx)u, Axu) = - I div(MCi;^) dx = 0 

2 Jr2 

for a smooth u G divu = 0. The general case follows by a 

standard approximation procedure. 

The next corollary gives the dissipative estimate for the solutions of 
fll.211 in TL^, and its proof is similar to that of fl2.8p . 
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Corollary 2.8. The strong solution satisfies the following estimate: 

< ||V.woire-“* + {2au)-%f, 

r^+T (2.18) 

V \\AMs)fds<u-^T\\gf + \\VMtW- 

In the following corollary we establish the smoothing property for 
the weak solutions of equation fll.2p . 

Corollary 2.9. Fort > 0, u{t) G and the following estimate holds: 

ll«(<)lll,. < cw (||«(0)||^ + llsf). 0 < « < 1, (2.19) 

where C = C{a,u). In particular, the attractor A is bounded in . 

Proof. Since ||Va:M(t)|P is integrable, for every 0 < T < 1, there exists 
T < T such that u{t) G FL^ and 

II V.«(r)f < 1 r ||V,«(s)f * < ^ (||«(0)||^ + a-T||9lO , 

where the second inequality follows from (12.81) . By (12.181) with the 
initial time f = r, we have 

l|V.«(r)f < ||v,«(T)f + (2aO-‘ll9ll" < Cr-‘ (||«(0)||^ + ||gf) . 
This and (12.8p prove (I2.19p and guarantee that the "HCball 
Bi := {no e H\ llnoll^i < R} 

is a bounded absorbing set for the solution semigroup S{t) if R is large 
enough. Thus, .4, C is bounded in FR and the proof is complete. □ 

Remark 2.10. Arguing analogously, it is not difficult to show that the 
smoothness of the global attractor A is restricted by the smoothness of 
the external forces g. In particular, the attractor will be C'°°-smooth 
(analytic) if the external forces are C°°-smooth (analytic). 

3. Fractal dimension of the attractor: upper bounds 

In this section, we prove that the global attractor constructed above 
has hnite fractal (box-counting) dimension and give explicit upper 
bounds for this dimension in terms of the physical parameters u and a 
and the norms of the right-hand side g in homogeneous Sobolev spaces. 
To this end, we hrst need to remind some dehnitions, see [37111] for more 
detailed exposition. 

Definition 3.1. Let A d FL be a compact set in a metric space FL. 
Then, by the Hausdorff criterion, for every e > 0 it can be covered by 
the hnite number of £-balls in R. Let Ae(4., FL) be the minimal number 
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of such balls. Then, the fractal (box-counting) dimension of ^ in "H is 
dehned via the following expression: 

In - 

The fractal dimension coincides with the usual dimension if the set A 
is regular enough (for instance, for the case where ^ is a Lipschitz 
manifold in "H), but may be non-integer for irregular sets (for instance, 
for the standard ternary Cantor set K C [0,1] this dimension is |^), 
see, for instance, [33] for more details. 

We now estimate the dimension of the global attractor A by means 
of the so-called volume contraction method PI, ip. The solution 
semigroup St is uniformly quasi-differentiable on the attractor in the 
sense that there exists (for a hxed f > 0) a linear bounded operator 
DS(t,uo) such that 

||S'(t)Mi - S{t)uo - DS{t,Uo) ■ {ui - Mo)|| < h{\\ui - MqII), (3.1) 

where h{r)/r —)■ 0 as r —)■ 0, and Uq, Ui G A. 

The quasi-differential DS{t,uo) is the solution operator ^ —)■ v{t) of 
the following equation of variations: 

dtv = L{t,uo)v := -Il[{v,Vx)u{t) + {u{t),Vx)v)-av+iyAa;V, u(0) = 

(3.2) 

where fl : [L^(M^)]^ —)■ "H is the Leray ortho-projection onto the diver¬ 
gence free vector helds, and u{t) = S{t)uo is the solution lying on the 
attractor and parameterized by uq G A. For the proof see [3|, where 
it is also shown that the solution semigroup is even differentiable for 
all Uq and the differential DS{t,uo) depends continuously on the 

point Uq. 

We dehne for m = 1,2... the numbers q{m) (the sums of the hrst 
m global Lyapunov exponents) 

^ pt rn 

g(m) = limsup sup sup - / y ^[L{t, Uq)v j, vj)dr, 

t^oo uo&A ^ Jo 

where the supremum closest to the integral is taken with respect to all 
L^-orthonormal families G . 

To dehne q{m) for all real m > 1 we just linearly interpolate between 
q{m) and q{m + 1) so that q{m) is now a piece-wise linear continuous 
function of m. 

The following theorem is the key technical tool for estimating the 
dimension of the attractor via the so-called volume contraction method. 
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Theorem 3.2. Let S(t) : "H —?• "H be the solution semigroup associated 
with problem fll. 2 p in a Hilbert space H and let A be a compact invariant 
set of St in H: S(T)A = A. 

Suppose that the semigroup S(t) is uniformly guasi-differentiable for 
every fixed t on A in the sense of fl3.1l) . 

Suppose further that the guasi-differential DS(t,uo) depends contin¬ 
uously on the initial point uo E A as a map DS(t, •) : mq —)■ 

Suppose that there exists number m > 0 such that q{m) < 0. Then 


dim/ ^ < m. 


For the proof of this theorem see |T3], |3^ in the case of Hausdorff 
dimension and [ 8 ] for the fractal dimension. 

Remark 3.3. The condition on the continuity of the quasi-differentials 
with respect to the initial point is redundant in the case of the Hausdorff 
dimension. It is also redundant in the case of the fractal dimension 
if the graph of q{m) lies below the straight line joining the points 
(m — l,g(m — 1 )) and (m,g(m)), where q{m) < 0 and q{m — 1 ) > 0 , 
see [To], E] . 

Also, in applications to inhnite dimensional dissipative dynamical 
systems an upper bound for q{m) is usually found in the form 


q{m) < —cim'^ + C 2 , 7 > 1. 


For example, as we shall shortly see, 7 = 1 in our case. In this case, 
also without the continuity condition, we have 


dim/ A < {c 2 /ciY^'^. 


To apply this theorem for obtaining the fractal dimension of the 
attractor A of the semigroup S{t) generated by damped Navier-Stokes 
equation fll.2p we need to state the Lieb-Thirring inequality which plays 
a fundamental role in estimating the quantities q{m). 

Lemma 3.4 (Lieb—Thirring inequality). Let {vjfjfi G be 

a family of orthonormal vector-functions and let divu/ = 0. Then the 
following inequality holds for p{x) = 



Proof. The proof (see [23]) is a reduction to the scalar case (which 
works in two dimensions) [ 8 ] and the use of the main result in [18] . □ 

We hrst set u = 1 and a = 1 and begin with estimating the m- 
trace of the operator L{t,uo) in the equation of variations fl3.2p for 
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this particular case. The general case will be reduced later to the case 
v = a = Ihy the proper scaling. 

Proposition 3.5. Let z/ = 1 and a = 1. Then, the following estimate 
holds: 


f*T ^ 


limsup sup sup 


T^-oo uoeA{vj}f^j^£H^T Jq “ 

1 


Uo)vj, Vj) dt < 


< —m + 


1 „ /TT limsup sup - / ||V^M(t)||^dt, (3.4) 

l0\/o T —^oo Jq 


where u{t) = S{t)uo. 


Proof. Let G P} be an orthonormal family in Pi. Then, 

integrating by parts and using that the vector fields Vi are divergence 
free, we get 


'^{L{u{t))vj, 


Vj) = 


i=i 


m m p 

j=i j=i -'R 


m 2 

j=l i,k=l 


m 2 


lit n 111 ^ 

-m-j2 IIv.uiir - 

i=l “^R^ j=l i,k=l 


x^u^Vj dx. 


To estimate the right-hand side above, we use the following point-wise 
inequality 

2 








(3.5) 


kA=l 


which holds for any v = (u^,u^) G and any Jacobi matrix = 
idxiU^)‘^i G such that + dx 2 u'^ = 0, see [H Lemma 4.1]. 

Indeed, setting 


v = {^,v), V^u = A-.= 


a b 
c —a 


we have by the Cauchy-Schwartz inequality 

{Av, v) = a{f‘^ - yf) + {b + c)fr] = a{f‘^ - yf) + ((6 -f c)/2) 2^t] < 
<\/of + {b + c)2/4 \/— rf'Y + = \/a? + {b + c)2/4 • |u|^ < 

<\/a? + {b"^ + c^)/2 ■ \v\^ = -y^/a? Pb"^ + + a? ■ |u|^, 

v2 
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and estimate (I3.5p is proved. Using this pointwise estimate and the 
Lieb-Thirring ineqnality for divergence free vector helds fl3.3p . we £- 
nally have 


m ^ 1 /* 

^{L{u{t))vj,Vj) <-m-^\\V:cVi\\^ + / P(a^)|V. 
i=i i=i V 2 Jm? 


u{x) I dx < 


< 


+ ^||p||||V^n(t)|| < 


2=1 


< 




< — m 


2=1 
ClT I 


%/2 


Vxu(t)\\ =-m + 


clt 11^ 
i=i 
1 


xVj\ 


1/2 




16\/3 




Integrating this ineqnality for t G [0, T] and taking the snpremnm over 
all Mo £ *^5 we obtain (13.dh and finish the proof of the proposition. □ 


Theorem 3.6. Let the assumptions of Theorem \2.1\ hold and let v = 
a = 1. Then the fractal dimension of the global attractor A of problem 
o satisfies the following estimate: 


where u(t) = S(t)uo. 

Proof. This estimate is a corollary of Theorem 13.21 and estimate fl3.4p . 

□ 


Thns, we only need to estimate the integral in the RHS of 03.61) . We 
assnme below that the right-hand side g belongs to the homogeneons 
Sobolev space = (—s G M with norm: 


\u\ 






(3.7) 


where 


HO ■= — f 

^ ^ 277 A 

is the Fonrier transform of u, see 


u{x)e dx 

for more details. Then, obviously, 
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and the following interpolation ineqnalities immediately follow from 
definition fl3.7p and the Canchy-Schwartz ineqnality: 

||w||^s < ||m|| llV^nll , 

—s II II — ll^aj^ll ll^3:;^ll 7 ^ ^ [^7 ^]* 

Corollary 3.7. Let the assumptions of Theorem, \2.1\ hold and let, in 
addition, g G H~^ for some s G [0,1] and n = a = 1. Then the fractal 
dimension of the attractor A of problem fll.2l) satisfies the following 
estimate: 

Proof. From the energy eqnality fl2.3p with u = a = 1, interpolation 
ineqnalities fl3.9p and the Yonng ineqnality, we get 


1 d 

2 dt 


I'^IIl 2 + 11^11^2 + ||Va;n||^2 — {g,u) < 


h-'*ll'^llh“ — 


— iiyiir/-» 


I 1 — 5 


< 


SP. 


I|Vxm||12||m||^ 2 —||y||^_. -r —1| Va7«||L2 


2 , Tuu2 


H-l|Va:W||L2 + —\\U\\l^, 

p q 


where p = ^, q = and positive nnmbers s, 6, 7 are snch that eSj = 1. 


2 2 
P<i = - 

Fixing now the parameter 7 in snch way that 

f 

q 

and integrating over f G [0,T], we arrive at 


r = i 


7 = 


_ 1 — 3 

1 — 2 


2T 


«(i)lt-M 0 )f)+ 1 - 


S’‘\ 1 


p J T 


|V.rM 


dt < 


From the dissipative estimate fl2.8p we conclnde that 


2 

H-=' 


lim snp — (||M(t)||^ - |1 m( 0)||2) = 0 
21 


and, therefore. 


lim snp—snp / ||Va;M 

T —^00 -1 JQ 


'dt < — {!- — 

p 


-1 


2 

H- = 


(3.11) 


and we only need to optimize the coefficient in the RHS with respect 
to £ and 5. Indeed, since = 1, we conclnde that 


£5 = 



e = 5 ^ 



1 —s 
2 
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and 


p 


-1 


1 /1-s 

2 


1-s 


a: (1 — |xs 


where x := 5^. Thus, it only remains to maximize the function 

/(x) := X (^1 - |x») 

on the interval x > 0: 


r// \ ^ S + 1 1 

/ (x) = 1- x» = 0 ^ X = 


s + l 


fix) = 


s + l V s + 1 


Inserting the obtained estimates into the right-hand side of (I3.1ip . we 
hnally get 


lim sup — sup / 

T —>-00 J 0 


I 


’ dt < 


1 + s 
1 — s 


2 

H-^ 


This estimate together with fl3.6p completes the proof of the corollary. 

□ 

The next corollary gives the analogous result for g G with s > 0. 


Corollary 3.8. Let the assumptions of Theorem \2.1\ hold and let, in 
addition, g E for some s G [0, 1] and u = a = 1. Then the fractal 
dimension of the attractor A of problem fll.2p satisfies the following 
estimate: 


dim f A < 


1 


1 + s 


(3.12) 


Proof. From the second energy equality (I2.17p with u = a = 1, inter¬ 
polation inequalities fl3.9p and the Young inequality, we get 


N^ll + II + II Aa;n|p = {g, A^u) < 
2 dt 

,|lV.n|n|A,n||i-* 


< 


^ \\y\\H 


< 


£2 

< — 
- 2 




+ s II A3;'U|| ^_s + 


+ yyllVxMlI^ H—^||Aa;np 


p q 

where the exponents p, q and the constants e, 6, 7 are exactly the same as 
in the proof of the previous corollary. Thus, using the strong dissipative 
estimate fl2.18p and arguing exactly as in the proof of the previous 
corollary, we end up with 


limsup—sup / WV xu{t)\f dt < 

T^oo T uo&A Jo 4 


1 — f 1 + s 


1 — s 
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which together with 


finishes the proof of the corollary. 


□ 


Thus, combining estimates fl3.10p and fl3.12|) . we end up with 


dimj A < 


1 — /I + Is 


64\/3 Vl - k 


hi 


(3.13) 


which holds for the case i> = a = lifgE H~^ fl and s G [—1,1]. 

Finally, we need the analogue of estimate fl3.13p for general z/, a > 0. 
We reduce this general case to the particular case i/ = a = 1 hj the 
proper scaling of t, x and u. Indeed, let u = u{t,x) be a solution of 
(I1.211 with arbitrary z/, a > 0. Then, taking 

1 


, , /a\i/2 

t := at, X = j X, 


u = 


tU 


{avYl‘^ 

we see that the function u{t',x') := -i-A^u{t,x) solves equation (11.21) 

with V = a = 1 and the external forces g{x') = a~^{av)~^^'^g{x). 
Since the fractal dimension of the attractor does not change under this 
scaling, using the obvious scaling properties of the i7®-norm: 


|2 

\h‘‘ 


\\9h-) 

we have proved the following result. 




2 


Theorem 3.9. Let the assumptions of Theorem \2.1\ hold (now with 
arbitrary positive v and a) and let, in addition g G FT® for some s G 
[—1,1]. Then, the fractal dimension of the attractor A in TL satisfies 
the following estimate: 


dim/ ^ < 


1 + |s| 


Id 


Z/N s 

a 


64\/3 o^z^^ '^aJ ^ ^ 

Proof. Indeed, the above estimate follows from fl3.13p and the identity 


ct^z/^ 


Z/N s 

a 


□ 


\i g & H ^ n H^, then the rate of growth of the estimate fl3.14l) with 
respect to z/ as z/ —)• 0 is the smallest when s = 1. In this case we have 

Corollary 3.10. Suppose that g G H^. Then the fractal dimension of 
the attractor A satisfies 

curl^fp 


1 

16\/3’ 




dim f A < 


(3.15) 


























21 


DAMPED NAVIER-STOKES EQUATIONS IN 

Proof. Since div^f = 0, it follows that 

ll^llii = l|Vx£/|P = II div^ll^ + II curlc/f = || curlc/f. 

□ 
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